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DESIGNING  OF  A  NONLINEAR  OPTIMAL  TERMINAL  GUIDANCE  LAW 

FOR  SPACE  INTERCEPTION 


Abstract:  This  article  concerns  a  nonlinear  kinematics  model  of  space  interception. 
Nonlinear  optimal  control  is  sought  through  use  of  inverse  systematic  methods.  A 
theoretical  nonlinear  terminal  guidance  law  of  global  state  feedback  is  advanced.  In 
addition,  a  simulation  test  is  carried  out  and  a  comparison  is  made  with  proportional 
guidance.  Results  show  that  this  guidance  law  performs  satisfactorily. 

Key  words:  nonlinear  guidance  law,  inverse  system  methods,  nonlinear  optimal  control, 
pseudolinear  system. 


GRAPHICS  DISCLAIMER 


All  figures,  graphics,  tables,  equations,  etc.  merged  Into  this 
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I. 


FORWARD 


Conventional  proportional  guidance  is  simple  and  easy  to 
carry  out.  However,  precision  is  relatively  low.  Moreover, 
line  of  sight  angular  velocity  divergence  phenomena  [4] 
exist.  For  this  reason,  there  is  a  need  to  search  out 
various  types  of  optimal  guidance  laws.  Due  to  actual 
interception  kinematics  models  being  nonlinear,  as  a  result, 
one  has  only  to  solve  optimal  control  problems  associated 
with  nonlinear  systems  and  only  then  is  it  possible  to 
obtain  relatively  high  precision  terminal  guidance  laws. 

With  regard  to  general  nonlinear  optimal  control  problems, 
people  still  have  not  found  a  unified  method  of  solution. 
However,  we  use  inverse  system  methods,  and  it  is  possible 
to  solve  this  concrete  nonlinear  optimal  control  problem 
encountered  in  terminal  guidance. 

II.  NONLINEAR  OPTIMAL  CONTROL  INVERSE  SYSTEM  METHODS 

Assume  the  general  form  of  nonlinear  systems  to 
be 


(2-1-a) 


(X=F{X,U)  X{t>)^X, 
*  W=G(X,U) 


(2-1-0) 

(2-1-6) 


(2-1-b) 

In  this,  U  is  an  m  dimension  input  vector.  Y  is  an  r 
dimension  output  vector.  X  is  an  n  dimension  state  vector. 


*  Numbers  in  margins  indicate  foreign  pagination. 
Commas  in  numbers  indicate  decimals. 


The  output  equation  (2-1-b)  is  also  capable  of  being  written 
as 


(yi  =  gi(X,U) 


(2-2) 


y,^gr{X,U) 


First,  carry  out  derivation  and  transformations  of 
equation  (2-2),  that  is,  in  sequence,  solve  for  the  relevant 
al,a2,...ar  order  derivatives  associated  with  time  period  t 
for  yl,y2,...yr,  constructing  the  new  equation 
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.yr''  =  hAX,U) 

hr^‘r^  =  hr{X,U) 


In  this,  numerical  values  for  ai(l^i^r)  are  defined  in 
accordance  with  the  forms  below 


(2-4-a) 


(2-4-b) 


,  eu 


k^Gi 


2 


In  this,  Fgt^(^^^gi^F 
,  and,  when  k  =  0,  one  has 


F^gi^F{F'‘-'gij 

F‘gi^gie 


Assuming  that,  during  the  process  of  the  derivatives 
described  above,  with  regard  to  l^i^r,  one  has  ai  <  <»,  in 
the  equation  set  (2-3),  if  U  can  act  as  a  function  of  x  and 
then,  it  is  possible  to  solve  for  the  explicit 
representation  as 


(2-5)  U=H-'(X,Y‘) 

In  this,  — .  Then,  it  is 

conveniently  possible,  from  equation  (2-1),  to  construct  a 
new  system  equation.  The  state  equations  can  be  obtained 
from  substituting  equation  (2-5)  into  equation  (2-1-a) .  The 
output  equation  is  then  form  (2-5).  As  a  result,  system 
/i;yu  equations  can  be  expressed  as 


(2-6-a) 


J-t  =  FCA%«-(A-,K')D,  A(/,)=A, 

(2-6-b) 
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It  can  be  demonstrated  [8]  that  the  system  expressed  by 
the  equations  described  above  is  the  inverse  systems  to 
system  (2-1). 

With  regard  to  the  initially  given  system  expressed  by 
equation  (2-1),  take  Ya  in  the  inverse  system  (2-6)  and 
replace  it  to  be  =  r-j.  -TO'  .  It  is  then 

conveniently  possible  to  obtain  the  a  order  integral  inverse 
system. 


(2-7-a) 

1  A'  =^FlX,H-'{X,<pn  A-('o)  =.V. 


(2-7-b) 


After  that,  before  taking  the  a  order  integral  inverse 
system  obtained  and  connecting  it  with  the  original  system, 
combined  systems  which  are  further  obtained  are  then 
corresponding  pseudo  linear  systems.  The  transmission 
relationship  can  be  expressed  as 


(2-8) 


diag{D^>,D‘‘‘,-,D'’')-Y  =  (p 


In  this,  Dai  is  the  ai  order  operator,  i  =  1,2, ...,r. 
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Pseudo  linear  systems  can  also  be  expressed  as  state  special 
forms 


ii—AZ  +  Bg) 

u=cz 


(2-9-b) 


Assuming  that  optimal  performance  indicators 
corresponding  to  system  (2-1)  are 


+±\'\xm)Q'X{t)+U^i)S'Umdt 

2  -  I, 


(2-10) 


Choosing  pseudo  linear  system  (2-9) 's  optimal 
performance  index  to  be 


(2-11) 


J{Z,(p)==Y^^(t,)RZ{t,) 

1 

+y],^  cz'(t)QZ(t)+<p^(t)s<p(t):idt 


If,  at  the  same  time  that  3{Z,(p)  min,  one  also  has 
J' (X,U)  min',  then,  it  is  possible  to  first,  with  regard  to 
pseudo  linear  (2-9)  and  in  accordance  with  solution  methods 
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for  linear  second  order  type  optimal  regulator  devices,  get 
pseudo  linear  system  optimal  control 
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(2-12) 


Assuming  the  relationship  between  state  vectors  X  and  Z 


(2-13) 


Z^TiX) 


comprehensive  consideration  of  forms  (2-7-b),  (2-12), 
as  well  as  (2-13),  it  is  finally  possible  to  obtain 
nonlinear  system  (2-1)  optimum  controls  which  are 


(2-14) 


=H-XX*,<p*iT(X*)  ,/)3 


III.  THE  SETTING  UP  OF  INTERCEPTION  DEVICE  AND  TARGET 
RELATIVE  MOTION  EQUATIONS 


Assuming  that  targets  move  along  elliptical  orbits,  in 
studying  instants  with  targets  at  point  M  and  interception 
devices  at  point  D,  the  relative  distance  between  the  two  is 
p.  Introducing  line  of  sight  coordinate  system  Dli7S»  , 

the  relationship  between  it  and  the  interception  device 
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three  dimensional  coordinate  system  Dxlylzl  is  determined  by 
altitude  angles  0  and  azimuth  angles  i|r  as  shown  in  Fig.l. 

The  mathematical  relationship  between  the  two  coordinate 
systems  is 


cosOcos^  sinfl  —cos6sinip  ' 

—  sin0cos^  cosd  sinOsiTpp 

Vi 

sintp  0  cos^ 

.  . 

Assuming  that  angular  velocities  of  rotation  of  line  of 
sight  coordinate  systems  with  respect  to  inertial  coordinate 
systems  are: 


^2-2)  S^Sj+S^+Sc 

In  this,  i3{,3fi.S;  are  projections  of  absolute  angular 
velocities  5  the  three  line  of  sight  coordinate  system 
axes.  When  target  orbit  altitudes  are  between  300km 
2500km,  considerations  of  terminal  guidance  laws  can  ignore 
differences  in  gravity.  Then,  in  studying  three  axis 
stability  plans,  and,  assuming  conditions  where  attitude 
controls  guarantee 
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(3-3) 


<ot=o 


relative  motion  equations  for  targets  and  interception 
devices  can  be  written  as : 


(3-4-a) 

,p  — pCcoS  +  oH  =0£ 

(3-4-b)  <2pc»£  +  P«>t  =  09 


(3-4-c) 

In  these,  .a{,ai7,a{  are  thrust  acceleration  a 

projections  on  the  three  line  of  sight  coordinate  system 
axes .  Due  to  the  forms  above  implicitly  containing  target 
movements,  this  article,  therefore,  hypothesizes  that 
targets  are  a  type  of  special  case.  Assume  that,  on 
interception  devices,  there  are  installed  four  normal  thrust 
engines  along  the  periphery  of  the  center  of  mass  as  shown 
in  Fig. 2.  That  is,  3,^=5  .  From  Fig.l,  it 

can  be  known  that: 
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Assume  that  attitude  control  systems  guarantee  0-*o,  ^>-*-0, 
Then,  one  has 


(3-6-a) 


(3-6-b) 


(3-6-c) 


,of=0 

ja,=o», 

'ot  =  a,, 


Moreover,  equation  (3-4)  changes  to  become 


(3-7-a) 

p(£ai  +  o|)  =0 
•j2pot  +  P“>C  =  0y, 

(3— 7— b)  Zpoii  pthii—Qg^ 


(3-7-c) 


The  equations  above  are  nothing  else  than  target  and 
interception  device  relative  motion  equations  in  three 
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dimensional  space.  Attitude  control  system  engine  fuel 
combustion  will  give  rise  to  changes  in  interception  device 
centers  of  mass.  Moreover,  once  four  normal  thrust  engines 
are  installed,  the  locations  are  not  easy  to  change. 
Therefore,  in  the  process  of  interception,  interception 
device  centers  of  mass  will  deviate  by  small  amounts  the 
planes  determined  by  the  centers  of  mass  of  the  four  normal 
direction  thrust  motors .  Assuming  that  the  errors  given 
rise  to  from  this  factor  as  well  as  the  errors  given  rise  to 
by  coupling  of  attitude  and  center  of  mass  movements  can  be 
ignored  in  engineering  terms,  then  kinematics  model  (3-7)  is 
reasonable.  From  this,  it  can  be  seen  that  model  (3-7)  is 
one  type  of  approximation  model.  This  type  of  approximation 
causes  the  model  to  be  very  greatly  simplified.  As  a 
result,  the  carrying  out  of  guidance  law  design  is 
facilitated. 

If  one  lets  X\=p,  Xi  =  p,  x,  —  o)n,  x,  —  coi>  U2  =  a>,» 

then,  equation  (3-7)  can  be  written  as  the  state 
equations 


(3-8-a) 

(3-8-b) 

(3-8-c) 

(3-8-d) 


.Vi(x5txJ) 

_  2^2X3 _ 

Xi  Xi 


2xtXi  Ut 

Xi  X, 


In  these,  xl  is  the  relative  distance 
between  interception  devices  and  targets.  x2  is  their 
relative  velocity.  x3  is  their  longitudinal  planar 
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angular  line  of  sight  velocity.  x4  is  transverse  planar 
angular  line  of  sight  velocity,  ul  and  u2  are, 
respectively,  accelerations  produced  in  association  with 
interception  device  thrusts  on  axes  Z1  and  Yl. 


IV.  NONLINEAR  OPTIMUM  GUIDANCE  LAW  DESIGN 


For  the  sake  of  convenience  in  the  study  of  guidance 
laws,  it  is  possible  to  take  equation  (3-8)  and  divide  it 
into  two  dimensional  models,  that  is,  longitudinal  planar 
relative  motion  equations 


(4-1-a) 


(4-1-b) 


.  _  2xtX,  .  til 
’  Xi  Xi 


(4-1-C) 

as  well  as  transverse  planar  relative  motion  equations 
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(4-2-a) 


(4-2-b) 


Xi  =  Xi 

it=x,x] 

■  _  2xtX4  .  tit 

X4  —  *T  ' 

Xi  Xi 


(4-2-c) 


Now,  taking  the  transverse  planar  case  in  order  to 
design  guidance  laws,  assume  system  (4-1) 's  output  equation 
is 
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(4-3) 


With  regard  to  equation  (4-3),  carry  out  derivation  and 
transformation  to  obtain 


.  .  2x,x,  u, 

(4-4)  ’  Jfi  X. 

After  transformation,  one  has: 
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(4-5) 


«i=  -Xx^i-2xtVt 


Let  V—ixi  ,  then,  pseudo  linear  system  numerical 
models  are: 


(4-6) 


are 


ix  =  <P 

Assume  that  pseudo  linear  system  state  variables 


(4-7) 

Then,  pseudo  linear  system  (4-6)  can  also  be  written  as 
state  controls  forms: 


(4-8-a) 

f^  =  9> 

hx=2 


(4-8-b) 
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If  one  takes  equation  (4-6)  and  substitutes  into 
equation  (4-5),  then,  one  has: 


(4-9) 


“i— 


After  considering  equations  (4-3)  and  (4-7),  equation 
(4,9)  is  changed  to  be: 


(4-10)  ti,— —x,q>-2x}Z 

From  equation  (4-10),  it  is  possible  to  write  pseudo 
linear  system  characteristic  indices  as 


(4-11 )  J Cgz’(0  -rsqi^(i)')dt 

In  the  equation,  r,  q,  and  s  are  weighting 
coefficients,  to  it  terminal  guidance  origination  time,  tb 
is  time  of  entry  into  guidance  blind  area.  From  equation 
(4-10)  and  equation  (4-11),  it  is  possible  to  see  that  one 
has  only  to  appropriately  select  weighting  coefficients. 
Then,  when  equation  (4-11)  reaches  a  minimum,  it  is  possible 
to  make  the  absolute  values  of  line  of  sight  angular 


16 


velocity  x3  and  control  ul  diminish.  That  is,  it  is 
possible  to  make  equation  (2-11),  which  corresponds  to  a 
matrix  associated  with  a  certain  set  of  weighting 
coefficients  R',Q',S',  reach  minimums.  As  a  result,  making 
connections,  it  is  possible  to  carry  on  in  accordance  with 
methods  described  in  Section  II. 


(4-12) 


In  this,  p(t)  satisfies  Riccati  scalar  quantity 
integral  equations  [ 7 ] : 


(4-13-a) 

fp(0=V(/)-9 

hiU)=r 

(4-13-b) 

Through  solving  Riccati  equations,  it  is  finally 
possible  to  obtain  pseudo  linear  system  optimum  control  as; 


(4-14) 


(r  +  Vsg)  +  (r  — *>/5g)exp  ^2 

- - - - -1± - -X*  if) . 

(r  +  's/sg)-(r-V5g)eip  J 
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Take  equation  (4-14)  and  substitute  into  equation 
(4-9).  One  then  obtains  the  longitudinal  planar  nonlinear 
optimum  guidance  law: 


(4-15) 


_  (r+ Vsg)  +  (r  —  V5g)exp  [2 

a?=  /— " - i= - x*Xt 

^  (r  +  VIg)  —  (r-'/sg)e'xp  1^2  /.^(f— ^j) j 


2x1x 


3  • 
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By  the  same  reasoning,  it  is  possible  to  obtain  the 
transverse  planar  nonlinear  optimal  guidance  law: 


(4-16)  1 


(r  +  >/5g)  +  (r  — N/sg)exp  [ 

- - - - - — + 

(r  +  */jg)-(r-Vag)exp  J 


Looking  at  the  two  equations  above,  this  type  of 
nonlinear  optimal  guidance  law  requires  the  carrying  out  of 
global  state  feedback  on  systems .  This  then  brings  with  it 
certain  difficulties  in  realizing  guidance  laws.  In 
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particular,  it  requires  the  measurement  of  distances. 
However,  we  are  capable  of  using  reconstructed  state  methods 
in  order  to  solve  this  problem.  In  summary,  although  the 
guidance  laws  we  put  forward  are  a  bit  complicated,  guidance 
performance  will,  however,  be  very  greatly  increased. 


V.  NUMERICAL  VALUE  SIMULATION 


Using  longitudinal  planar  nonlinear  guidance  laws  as  an 
example,  simulation  tests  were  carried  out.  When 

x,(fo)  =  100km, x,(fo)  = -lOkm/sec,  5i-,(fo)  =5 X  10  ^rad/sec  , 

simulation  results  for 

guidance  law  (4-15)  are  as  shown  in  Fig.  3,  4,  and  5  (among 
these,  selections  are  made  for  r=10,  q=60,  s=110,  to  =  0 
sec,  and  tb  =  10  sec).  Moreover,  proportional  guidance  law 
simulation  results  are  seen  in  Fig. 6  and  7.  In  accordance 
with  calculation  formulae  for  the  amounts  by  which  targets 
are  missed  [ 9 ] , 


(5-1) 


Pinlo 


Xtb^  xit+y,* 

V  xlb  +  xUixU  +  xU) 


it  is  possible  to  solve  for  the  amounts  by  which 
targets  are  missed  as  associated  with  nonlinear  guidance  law 
(4-15)  and  conventional  proportional  guidance  laws  to  be, 
respectively,  5.4x10-7  m  and  10.32  m.,  assuming  the  aiming 
blind  area  range  is  2000m.  In  nonlinear  guidance  laws,  the 
*  values  of  weighting  coefficients  r,  q,  and  s  influence 
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guidance  characteristics.  Through  simulations,  it  is 
possible  to  know  that,  the  larger  the  value  of  q/s  is,  the 
better  guidance  performance  is.  Moreover,  the  influence  of 
the  value  of  r  on  guidance  characteristics  is  very  tiny. 
Simulation  results  also  tell  us  that  the  absolute  values  of 
control  signals  associated  with  nonlinear  guidance  law 
(4-15)  are  monotonic  reduction  functions  of  time.  Moreover, 
they  will  not  exceed  maximum  thrust  ranges  associated  with 
interception  device  motors.  Besides  this,  conventional 
proportional  guidance  indeed  has  line  of  sight  angular 
velocity  divergence  phenomena.  Moreover,  as  far  as  the 
nonlinear  guidance  laws  discussed  above  are  concerned,  the 
line  of  sight  angular  velocities  drop  monotonically  and 
converge  to  zero.  In  summary,  the  performance  of  the 
nonlinear  guidance  laws  given  in  this  article  are  indeed 
very  good. 


1»XO 


Fig. 3  Changes  in  Relative  Distances  Under  the  Control 
of  Nonlinear  Guidance  Laws 


Fig. 4  Changes  in  Line  of  Sight  Angular  Velocities 
Under  Nonlinear  Guidance  Law  Control 


Fig. 5  Change  States  Associated  with  Control  Signals  in 
Nonlinear  Guidance  Laws 
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/24* 


T(Sec) 


Fig. 6  Changes  in  Relative  Distance  Under  Proportional 
Guidance  Law  Control 


Fig. 7  Changes  in  Line  of  Sight  Linear  Velocities 
Under  Proportional  Guidance  Law  Control 


Numbers  in  margins  indicate  foreign  pagination. 

Commas  in  numbers  indicate  decimals.  « 
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VI.  CONCLUSIONS 


♦ 


This  article  makes  use  of  inverse  system  methods  and 
studies  nonlinear  optimal  control  problems.  In  conjunction 
with  that,  it  puts  forward  nonlinear  optimal  guidance  laws. 
The  guidance  laws  in  question  furnish  a  theoretical 
foundation  for  increasing  space  intercept  terminal  guidance 
performance.  Of  course,  during  actual  realization,  there 
still  exist  a  number  of  problems — for  instance,  problems 
with  obtaining  relative  distance  signals  xl(t)  in  it. 
Please  see  reference  [6].  This  article  will  not  provide 
superfluous  details . 
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